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Abstract 

We summarize the key results of a recent global analysis of inclusive single- 
charged-hadron production in high-energy colliding-beam experiments. In the frame- 
work of the parton model of quantum chromodynamics at next-to-leading order 
(NLO), fragmentation functions (FFs) for charged pions, charged kaons, and (anti) 
protons are extracted from experimental data of e + e~ annihilation at the Z-boson 
resonance and at centre-of-mass energy = 29 GeV. This fit also yields a new 
NLO value for the strong-coupling constant, namely a4 '(Mz) = 0.1170 ± 0.0073. 
The scaling violations encoded in the FFs through the Altarelli-Parisi evolution 
equations are tested by confronting e + e~ -annihilation data from DESY DORIS, 
DESY PETRA, and CERN LEP2 with NLO predictions based on these FFs. Com- 
parisons of pp data from CERN SppS and the Fermilab Tevatron, 7p data from 
DESY HERA, and 77 data from LEP2 with the corresponding NLO predictions 
allow us to test the universality of the FFs predicted by the factorization theorem. 

PACS numbers: 13.60.Le, 13.65.+i, 13.85.Ni, 13.87.Fh 

1 Introduction 



In the framework of the QCD-improved parton model, the inclusive production of single 
hadrons is described by means of fragmentation functions (FFs), D^(x,fi 2 ). The value of 
D%(x, fi 2 ) corresponds to the probability for the parton a produced at short distance l//x to 
form a jet that includes the hadron h carrying the fraction x of the longitudinal momentum 
of a. Unfortunately, it is not yet possible to calculate the FFs from first principles, 
in particular for hadrons with masses smaller than or comparable to the asymptotic 
scale parameter, A. However, given their x dependence at some energy scale fi, the 
evolution with fi may be computed perturbatively in QCD using the timelike Altarelli- 
Parisi equations [p]]. This allows us to test QCD quantitatively within one experiment 
observing single hadrons at different values of centre-of-mass (CM) energy y/s (in the case 
of e + e~ annihilation) or transverse momentum pt (in the case of scattering). Moreover, 
the factorization theorem guarantees that the D^(x,fi 2 ) functions are independent of the 



process in which they have been determined, and represent a universal property of h. 
This enables us to make quantitative predictions for other types of experiments as well. 

After the pioneering leading-order (LO) analyses of pion, kaon, and charmed-meson 
FFs in the late 1970s |§, there had long been no progress on the theoretical side of 
this field. In the mid 1990s, next-to-leading-order (NLO) FF sets for 7T°, tt ± , K±, and 
7] mesons were constructed through fits to data of e + e~ annihilation, mostly generated 
with Monte Carlo event generators ||. In 1994/95, the author, in collaboration with 
Binnewies and Kramer, extracted 7r =1= and K ± FFs through fits to SLAC-PEP and partially 
preliminary CERN-LEP1 data and thus determined the strong-coupling constant to be 
af\M z ) = 0.118 (0.122) at NLO (LO) || (BKK). However, these analyses suffered from 
the lack of specific data on the fragmentation of tagged quarks and gluons to tt^^, K ± , 
and p/p hadrons. 

During the last five years, the experiments at LEP1 and SLAC SLC have provided 
us with a wealth of high-precision information on how partons fragment into low-mass 
charged hadrons, so as to cure this problem. The data partly comes as light-, o, and b- 
quark-enriched samples without j5|H[7| or with identified final-state hadrons (7r ± , K^, and 
p/p) [j§|,j9l,|10] or as gluon-tagged three-jet samples without hadron identification |fTT|J]~^JT3|1 . 
Motivated by this new situation, the author, together with Kramer and Potter, recently 
updated, refined, and extended the BKK |4j] analysis by generating new LO and NLO sets 
of 7r ± , K ± , and p/p FFs fl^.Q By also including in our fits 7r ± , K ± 1 and p/p data (without 
flavour separation) from PEP ||15|| , with CM energy = 29 GeV, we obtained a handle 
on the scaling violations in the FFs, which enabled us to determine the strong-coupling 
constant. We found af\M z ) = 0.1181 ± 0.0085 at LO and af\M z ) = 0.1170 ± 0.0073 
at NLO |16| |. These results are in perfect agreement with what the Particle Data Group 
currently quotes as the world average, af\M z ) = 0.1181 ± 0.002 |T7[. Our strategy was 
to only include in our fits LEP1 and SLC data with both flavour separation and hadron 
identification |]8|,|9|,p!0| , gluon-tagged three-jet samples with a fixed gluon-jet energy |TT|Jl2]| , 
and the n^, K^, and p/p data sets from the pre-LEPl/SLC era with the highest statistics 
and the finest binning in x |L5]. Other data served us for cross checks. In particular, we 
probed the scaling violations in the FFs through comparisons with 7r ± , K , and p/p data 
from DESY DORIS and DESY PETRA, with CM energies between 5.4 and 34 GeV @. 
Furthermore, we tested the gluon FF, which enters the unpolarized cross section only at 
NLO, by comparing our predictions for the longitudinal cross section, where it already 
enters at LO, with available data |]5|JT9|. Finally, we directly compared our gluon FF with 
the one recently measured by DELPHI in three-jet production with gluon identification 
as a function of x at various energy scales fi |13|]. All these comparisons led to rather 



encouraging results. We also verified that our FFs satisfy reasonably well the momentum 
sum rules, which we did not impose as constraints on our fits. 



Very recently, we extended our previous tests of scaling violations [|14J] to higher energy 
scales by confronting new data on inclusive charged-hadron production in e + e~ annihila- 



1 A FORTRAN subroutine which retur ns the values of the D^(x : n 2 ) functions for given values of x 
and [J? may be downloaded from the URL http://www.desy.de/~poetter/kkp.html or obtained upon 
request from the authors. 
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tion from LEP2 PO 



based on our FFs 21 



with ranging from 133 GeV up to 189 GeV, with NLO predictions 
Furthermore, we quantitatively checked the universality of our 



FFs by making comparisons with essentially all available high-statistics data on inclusive 
charged-hadron production in colliding-beam experiments [^TJ. This includes pp data 
from the UAl and UA2 Collaborations pj at CERN SppS and from the CDF Collabo- 



ration [p3] at the Fermilab Tevatron, 7p data from the HI and ZEUS Collaborations [2| 
at DESY HERA, and 77 data from the OPAL Collaboration |2| at LEP2. 
alternative sets of NLO FFs for ^ 



In 2000, alternative sets of NLO FFs for 7r ± , K ± [26], and charged hadrons [[26,27] have 
become available. They are based on different collections of experimental data and on 
additional theoretical assumptions. In Ref. ||26|| , power laws were assumed to implement 
a hierarchy among the valence- and sea-quark FFs. In Ref. |27| , the renormalization and 
factorization scales were identified and adjusted according to the principle of minimal 
In order to estimate the present systematic uncertainties in the FFs, we 



sensitivity [28 



compare NLO predictions for inclusive charged-hadron production consistently evaluated 
with the three new-generation FF sets |1^J26| , |27]| . 

In this contribution, we summarize the key results obtained in Refs. [n| , |TT)| , |2~T] | . In 
Section Q, we present some details of our global fits and assess the quality of the resulting 
FFs. In Section ||, we discuss the determination of af\Mz) from the scaling violations in 
the FFs and compare our NLO result with those of different determinations. In Section |4], 
we present comparisons of our NLO predictions for inclusive charged-hadron production 
with e + e~ data from LEP2 [|2(J , pp data from SppS []22j and the Tevatron [ 53J , 'jp data from 
HERA ||24j|, and 77 data from LEP2 fl25|. Our conclusions are summarized in Section |5L 



2 Determination of the FFs 

The NLO formalism for extracting FFs from measurements of the cross section da/dx 
of inclusive hadron production in e + e~ annihilation was comprehensively described in 
Ref. [|J. We work in the MS renormalization and factorization scheme and choose the 
renormalization scale /i and the factorization scale Mf to be ^ = Mf = ^y^, except 
for gluon-tagged three-jet events, where we put /1 = Mf = 2^E jet , with E )et being the 
gluon jet energy in the CM frame. Here, the dimensionless parameter £ is introduced 
to determine the theoretical uncertainty in af\Mz) from scale variations. As usual, 
we allow for variations between £ = 1/2 and 2 around the default value 1. For the 
actual fitting procedure, we use x bins in the interval 0.1 < x < 1 and integrate the 
theoretical functions over the bin widths as is done in the experimental analyses. The 
restriction at small x is introduced to exclude events in the region where mass effects and 
nonperturbative intrinsic-transverse-momentum effects are important and the underlying 
formalism is insufficient. On the other hand, our analysis should be rather insensitive to 
nonperturbative effects at x values close to 1, since the experimental errors are very large 
there. We parameterize the x dependence of the FFs at the starting scale fio as D%(x, /Xq) = 
Nx a (l— x) 13 . We treat N, a, and j3 as independent fit parameters. In addition, we take the 
asymptotic scale parameter A^L appropriate for five quark flavors, as a free parameter. 
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Table 1: CM energies, types of data, and Xdf values obtained at LO and NLO for the 
various data samples. 
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Thus, we have a total of 46 independent fit parameters. The quality of the fit is measured 
in terms of the x 2 value per degree of freedom, Xdfi for all selected data points. Using a 



multidimensional minimization algorithm [2S], we search this 46-dimensional parameter 



space for the point at which the deviation of the theoretical prediction from the data 
becomes minimal. 

The Xdf values achieved for the various data sets used in our LO and NLO fits may 
be seen from Table [I]. Most of the Xdf values lie around unity or below, indicating 
that the fitted FFs describe all data sets within their respective errors. In general, the 
Xbf values come out slightly in favor for the DELPHI data. The overall goodness of 
the NLO (LO) fit is given by Xbf = 0-98 (0.97). The goodness of our fit may also be 
judged from Figs. |TJ(a) and [Tj(b) , where our LO and NLO fit results are compared with 
the ALEPH §0, DELPHI @, OPAL JTJ, and SLD data. In Fig. |(a), we study 
the differential cross section (1/ 'a tot )da h / 'dx for n ± , K ± , p/p, and unidentified charged 
hadrons at y/s = 91.2 GeV, normalized to the total hadronic cross section 0"t o t, as a 
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function of the scaled momentum x = Iphj^fs. As in Refs. [§,[nj, we assume that the 
sum of the 7r ± , K , and p/p data exhaust the full charged-hadron data. We observe 
that, in all cases, the various data are mutually consistent with each other and are nicely 
described by the LO and NLO fits, which is also reflected in the relatively small Xdf values 
given in Table |1]. The LO and NLO fits are almost indistinguishable in those regions of 
x, where the data have small errors. At large x, where the statistical errors are large, 
the LO and NLO results sometimes moderately deviate from each other. In Fig. |l](b), we 

measurements of the gluon FF in gluon-tagged 
26.2 and 40.1 GeV, respectively, with our LO 



compare the ALEPH || and OPAL fl 
charged-hadron production, with E- Je t 
and NLO fit results. The data are nicely fitted, with Xdf values of order unity, as may 
be seen from Table p]. By the same token, this implies that these data sets | Tiyi2" 
mutually consistent. 



arc 
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(a) (b) 

Figure 1: (a) Comparison of data on inclusive charged-hadron production at -^s = 
91.2 GeV from ALEPH § (triangles), DELPHI (circles), and SLD |I| (squares) 
with our LO (dashed lines) and NLO (solid lines) fit results. The upmost, second, third, 
and lowest curves refer to charged hadrons, ir ± , K^, and p/p, respectively, (b) Compar- 
ison of three-jet data on the gluon FF from ALEPH |lTJ with E- Jet = 26.2 GeV (upper 
curves) and from OPAL Jl2[ with E^ et = 40.1 GeV (lower curves) with our LO (dashed 
lines) and NLO (solid lines) fit results. 
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3 Determination of ct^(Mz) 



Since we included in our fits high-quality data from two very different energies, namely 29 
and 91.2 GeV, we are sensitive to the running of ct s (/j,) and are, therefore, able to extract 
values of A^. We obtain = 88+^ + 2 3 3 MeV at LO and A§L = 213+^9 MeV at 
NLO, where the first errors are experimental and the second ones are theoretical. The 
experimental errors are determined by varying A^ in such a way that the total Xdf value 
is increased by one unit if all the other fit parameters are kept fixed, while the theoretical 
errors are obtained by repeating the LO and NLO fits for the scale choices £ = 1/2 and 
2. From the LO and NLO formulas for ai n/ ^(yu), we thus obtain 

a f(M z) =0.1181 +"-0"58 +0.0006 
s y ' -0.0069-0.0049 v ; 

a f(M z) =0.1170 +a0055+a0017 (NLO), (1) 
s K ' -0.0069-0.0025 ^ h K J 

respectively. As expected, the theoretical error is significantly reduced as we pass from LO 
to NLO. Adding the maximum experimental and theoretical deviations from the central 
values in quadrature, we find = (88 ± 41) MeV and af\M z ) = 0.1181 ± 0.0085 at 

LO and A^ = (213 ± 79) MeV and af\M z ) = 0.1170 ± 0.0073 at NLO. We observe 
that our LO and NLO values of af'{M z ) are quite consistent with each other, which 
indicates that our analysis is perturbatively stable. The fact that the respective values of 
A^ significantly differ is a well-known feature of the MS definition of ai n/ ^(/i) j30] . 
Our values of A^ and af\M z ) perfectly agree with those presently quoted by the 

Particle Data Group (PDG) [[TJ as world averages, AJJ = 208 +| MeV and af\hl z ) = 
0.1181 ± 0.002, respectively. Notice that, in contrast to our LO and NLO analyses, the 
PDG evaluates A^ from a^(M z ) using the three-loop relationship [J30| J. The PDG 
combines twelve different kinds of (M z ) measurements, including one from the scaling 
violations in the FFs |3l|, by minimizing the total \ 2 value and thus obtains af\M z ) = 
0.1181 ± 0.0014 with \ 2 — 3.8. The world average cited above is then estimated from the 
outcome by allowing for correlations between certain systematic errors. It is interesting 
to investigate how the world average of af\M z ) is affected by our analysis. If we replace 
the value af\M z ) = 0.125 ± 0.005 ± 0.008 resulting from previous FF analyses |3lfl , 
which enters the PDG average, with our new NLO value, then we obtain af\M z ) = 
0.1180 ± 0.0014 with x 2 — 3.22, i.e., the face value of the world average essentially goes 
unchanged, while the overall agreement is appreciably improved. This is also evident 
from the comparison of Fig. ^, which summarizes our updated world average, with the 
corresponding Fig. 9.1 in Ref. [17]. We observe that the central value of our new NLO 
result for af\M z ) falls into the shaded band, which indicates the error of the world 
average, while in Fig. 9.1 of Ref. [17 the corresponding central value |3l| exceeds the 
world average by 3.5 standard deviations of the latter, which is more than for all other 
eleven processes. Furthermore, our new NLO result has a somewhat smaller error (0.0073) 
than the corresponding result [31] used by the PDG (0.009). This is due to a marked 
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decrease in the theoretical error, which may be attributed to a different choice of input 
data, especially at low CM energies. If we take the point of view that our new NLO value 
of a^(Mz) should rather be combined with the result from the previous FF analyses |H] 
before taking the world average, then the latter turns out to be af\M z ) = 0.1181±0.0014 
with x 2 = 3.34. 
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Figure 2: Summary of the values of af\Mz) from various processes. The errors shown 
represent the total errors including theoretical uncertainties. 



4 Global Analysis of Collider Data 



We now review recent tests of the scaling violations and the universality of our FFs |21 



On the one hand, we confronted new data on inclusive charged-hadron production from 
LEP2 PDf , with CM energies ranging from 133 GeV to 189 GeV, with our NLO predictions. 
On the other hand, we performed a global NLO analysis of essentially all high-statistics 
data on inclusive charged-hadron production in colliding-beam experiments, including 
pp scattering at SppS |22| and the Tevatron f23|, scattering at HERA |24]], and 77 



scattering at LEP2 |25[. In the cases of hadroproduction and photoproduction, we set 
fi = Mf = ^pt- As for the parton density functions (PDFs) of the proton, we employ set 
CTEQ5M provided by the CTEQ Collaboration with A^ = 226 MeV. As for the 
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photon PDFs, we use the set by Aurenche, Fontannaz, and Guillet (AFG) |33|]. In the 
following, we always consider the sum of positively and negatively charged hadrons. 

In all cases, we found reasonable agreement between the experimental data and our 
NLO predictions as for both normalization and shape, as may be seen from Fig. |3|. The 
majority of the data sets are best described with the central scale choice £ = 1. Exceptions 
include the UAl data sets with \fs = 200 and 630 GeV and the UA2 data set with 
1 < \y\ < 1.8, which prefer £ = 1/2, as well as the ZEUS data, which favours £ = 2. 
However, if we estimate the theoretical uncertainty due to unknown corrections beyond 
NLO by varying £ between 1/2 and 2, as is usually done, then it is justified to state 
that all the considered data sets agree with our NLO predictions within their errors. We 
hence conclude that our global analysis of inclusive charged-hadron production provides 
evidence that both the predicted scaling violations and the universality of the FFs are 
realized in nature. 



5 Conclusions 



Owing to the high-statistics experimental information on how partons fragment into low- 
mass charged hadrons provided by the LEP1 and SLC experiments, the determination of 
NLO FFs advanced to a level of precision which is comparable to the one familiar from 
similar analyses for PDFs. In this presentation, we reviewed recent LO and NLO analyses 
of tt^, K^, and p/p FFs [fL4| , which also yielded new values for af\M z ) pf. Although 
these FFs are genuinely nonperturbative objects, they possess two important properties 
that follow from perturbative considerations within the QCD-improved parton model and 
are amenable to experimental tests, namely scaling violations and universality. The scaling 
violations were tested [[14],[n] by making comparisons with data of e + e~ annihilation at 
CM energies below and above |20] those pertaining to the data that entered the fits. 



The universality property was checked |2T| by performing a global study of high-energy 
data on hadroproduction in pp collisions p2|j23[| and on photoproduction in e ± p [0] and 
collisions. 



e T e 
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As is well known, the gluon FF enters the prediction for the unpolarized cross section 
of inclusive hadron production in e + e~ annihilation only at NLO, while at LO it only 
contributes indirectly via the /i 2 evolution. In order to nevertheless have a handle on it, 
we included in our fits jl4| experimental data on gluon-tagged three-jet events from LEP1 
11.12]. Furthermore, we checked that our predictions for the longitudinal cross section, 



where it already enters at LO, agree well with available data [f|,|T9|. On the other hand, 
the gluon FF is known to play a crucial role for pp, 7p, and 77 scattering at low values 
of px- Thus, the comparisons performed here provide another nontrivial test of the gluon 
FF. 

As we have seen in Fig. |3|, the theoretical uncertainty of the NLO predictions due to 
scale variations significantly decreases as px increases. In order to perform more meaning- 
ful comparisons, it would, therefore, be desirable if pp, jp, and 77 experiments extended 
their measurements out to larger values of pt- Furthermore, in order to render such com- 
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parisons more specific, it would be useful if these experiments provided us with separate 
data samples of 7r ± , K^, and p/p hadrons. 

We also estimated the current systematic uncertainties in the FFs by comparing our 
NLO predictions for flavour-tagged inclusive charged-hadron production with those ob- 
tained from two other up-to-date FF sets p6|p7| . Apart from a difference in the 6-quark 
FFs at medium to large values of x, which may be traced to the incompatibility of two 
underlying b- quark-specific data samples [0,|], all three FF sets |[T4| . [26| . ^7| mutually agree 
within the present experimental errors. 
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Figure 3: Comparisons of data on inclusive charged- hadron production in (a) e + e anni- 
hilation at LEP2 with y/s = 133, 161, 172, 183, and 189 GeV (from bottom to top in this 
order), (b) pp scattering at SppS, (c) jp scattering at HERA, and (d) 77 scattering at 
LEP2 integrated over 77-invariant-mass intervals 10 < W < 30 GeV, 30 < W < 55 GeV, 
55 < W < 125 GeV, and 10 < W < 125 GeV (from bottom to top in this order) with 
NLO predictions based on our FFs. 
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